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Abstract. Several problems of interest in planetary infrared remote sensing are 
investigated using a new radiative-conductive model of energy transfer in regoliths. the 
solid-state greenhouse effect, thermal beaming, and reststrahlen spectra. The results of the 
analysis are as follows: (1) The solid-state greenhouse effect is self-limiting to a rise of a 
few tens of degrees in bodies of the outer solar system. (2) Non-Lambertian directional 
emissivity can account for only about 20% of the observed thermal beaming factor. T e 
remainder must have another cause, presumably surface roughness effects. (3) The 
maximum in a reststrahlen emissivity spectrum does not occur exactly at the Christiansen 
wavelength where, by definition, the real part of the refractive index equals one, but 
rather at the first transition minimum in reflectance associated with the transition from 
particle scattering being dominated by volume scattering to that dominated by strong 
surface scattering. The transparency feature is at the second transition minimum and does 
not require the presence of a second band at longer wavelength for its occurance. 
Subsurface temperature gradients have only a small effect on emissivity bands. 


1. Introduction 

Hapke [this issue] (hereinafter refered to as paper 1) intro- 
duced a new model for thermal and radiative transfer in plan- 
etary regoliths. This model simultaneously solves the radiative 
transfer equation for incident visible (including the near-UV, 
visible, and near-IR wavelength region) and emitted thermal 
infrared radiation and the heat conduction equation. Approx- 
imate analytic expressions are obtained for the distribution of 
visible and IR radiation fluxes and the temperature as a func- 
tion of depth below the surface of the medium, the visible 
hemispherical reflectance, and IR directional and hemispher- 
ical emissivities of a medium of an iso tropically scattering par- 
ticles. The solid-state greenhouse effect and the contribution of 
radiation to the thermal conductivity appear intrinsically in the 
model with no ad hoc additions. 

In this paper the new model is applied to three problems of 
interest in the remote sensing of planetary regoliths. Only 
time-independent solutions are considered. Section 2 analyzes 
the solid-state greenhouse in outer planet satellites. Section 3 
treats the effect of a non-Lambertian emissivity on thermal 
beaming. Section 4 discusses the relation between the spectra 
of the refractive index and emissivity in a reststrahlen band. 

The notation used in this paper is identical to that in paper 

1, and all the equations are taken exactly from that reference, 
which should be consulted for notation and definitions. 

2. The Solid-State Greenhouse Effect 

An atmospheric greenhouse effect can occur on a planet 
whose atmosphere is transparent to visible light but which 
absorbs thermal IR radiation. The energy of visible sunlight is 
absorbed at the ground level but is radiated back into space as 
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IR radiation from upper layers of the atmosphere. This re- 
quires a temperature gradient to transport heat from the sur- 
face to the radiating layer. On Venus this effect produces a 
surface temperature that is nearly triple the radiative equilib- 
rium temperature. Brown and Matson [1987] and Matson and 
Brown [1989] have pointed out that a strong solid-state green- 
house effect could occur in a planetary regolith that is weakly 
absorbing in the visible and strongly absorbing in the IR and 
that this would have important implications for the interior 
temperatures of icy satellites. However, they described the 
visible radiation source term by approximate, empirical expres- 
sions, rather than using the radiative transfer equation, so that 
their results are quantitatively uncertain. 

The results of the radiative equilibrium calculation of paper 
1 will be used to calculate the subsurface temperature rise in a 
regolith in equilibrium with sunlight. Equation (58) of paper 1 
for the reduced temperature (temperature relative to the 
blackbody radiative equilibrium temperature T r ) at infinite 
optical depth below the surface is 


T*(°°) - T*~ t rqT*\™) 


+ q 


B 


tt<Pt.s + K [ Jv VovA + 4 


( 1 ) 


where all quantities are defined in paper 1; (p* Ts ,A, and B are 
given by equations (56), (50b), and (50c) of paper 1, respec- 
tively, T * is the reduced surface temperature, given by (61) and 
(56) of paper 1; q = J JE* r kT r = cr 0 T; 3 / nE* r k, J v is incident 
visible irradiance (including the near-UV, visible, and near- 
IR), cr 0 is Boltzmann’s constant, T r is blackbody radiative equi- 
librium temperature, E% is reduced extinction coefficient of 
the regolith, and k is solid-state thermal conductivity. 

It is instructive to first calculate the temperature for the 
relatively simple case where the particles of the regolith are 
large compared to the wavelength and scatter isotropically. 
Then several of the quantities that occur in the parameters in 
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(1) can be simplified by putting Q E „ = Q ET = 1, fi v = fi T = 
0, E„ = E r , and = 1. After substituting the 

appropriate expressions from paper 1 for .4, B , and <p^ v , the 
quantities inside the square brackets in (1) can be factored and 
combined to yield 

r*(°°) ^ T* ~ TrqT*\cc) + q ^(\ + 2/x 0 )/2, (2) 


where ju, 0 = cos i and i is angle of incidence. 

In both (1) and (2) the first term on the right-hand side is the 
reduced surface temperature, while the second term accounts 
for heat transfer by thermal radiation. The third term describes 
the relative rise in subsurface temperature, that is, the green- 
house effect. Note that for the simplified conditions of (2) this 
term is independent of both the visible and thermal constants 
of the medium and is directly proportional to the incident 
irradiance J v . It has its maximum value of 3g/2 when jx 0 — 1. 
At the distance of Jupiter, a regolith having a high single 
scattering albedo, but otherwise like lunar soil, would have q — 
0.023 (see paper 1), so that 3q/2 - 0.035, corresponding to 
a greenhouse effect of only a few percent. The reason the 
relative temperature rise is so small is that although a small 
absorption coefficient allows the visible radiant power to be 
deposited deep in the regolith, it also causes a high visual 
albedo, so that less power enters the medium. The two effects 
exactly cancel each other when the particles are isotropic scat- 
tered. 

Hence small absorption in the visible combined with high IR 
opacity is not, by itself, sufficient to cause a large solid-state 
greenhouse effect in a regolith. However, the temperature rise 
can be increased if the particles are forward scattering, since 
then both the absorbance and the albedo can be small in the 
visible. To calculate the rise in temperature in this case, the full 
equation (1) from (58) of paper 1 must be used. Substituting 
from paper 1, factoring, collecting terms, and reverting to 
unstarred (unreduced) quantities gives 
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angle of incidence at the time of the measurements was \x Q = 
900/2800 = 0.32. 

The relative flux of visible light inside the snowpack is de- 
scribed by equation (49) of paper 1 for ip * . This equation 
consists of two terms. The first term describes the incident 
unscattered sunlight and is negligible deeper than a few par- 
ticle diameters. Below this depth only the second term re- 
mains, which describes the multiply scattered sunlight and is 
proportional to exp (-2y*r*). Hence x z can be identified 
with 2y*r* = 2£ v y v Ejz, so that \ = 2£ T /y„is„. 

Now, E v ~ NcrQ Ev , where N is number of particles per unit 
volume, cr is particle cross sectional area, and Q Ev is extinction 
efficiency. Since the particles are large compared with both 
visible and thermal radiation, Q Ev — Q E t ~ 1- The bulk 
density is p = p s .NTrd 3 /6, so N = 6<j>/iTd 3 , where </> = p/p s 
is the filling factor; since p, — 1, </> — 0.30. Thus E v = 
A 7 ctQ E j , “ (64>/Trd 3 ) * (m/ 2 /4) * 1 = 3<#>/2 d, so that xd = 
2£ v y v (3<i>/2d)d = 0.03, and f w y„ = 0.033. This gives one 
relation between 0„ and w 7/ . A second relation involving 0 T , 
and h»„ is the measured albedo of the snow, A h = 0.83 at 
/x, j = 0.32, which is given by equation (65) of paper 1. The two 
expressions can be solved simultaneously to give w v - 0.997 
and 0„ = 0.672, so that £„ — 0.574 and y v = 0.0574. 

The IR optical properties of snow are poorly known. Irvine 
and Pollack [1968] estimate on the basis of highly uncertain 
data that in the 10-25 pm wavelength range a perfect sphere 
of water ice with radius 10 pm has w — 0.8 and has an average 
cosine of the scattering angle —0.7. However, Salisbury et al. 
[1994a] measured the IR reflectance of water frost to be less 
than 1%, which implies that w must be very small there. Thus 
y r cannot be greatly different from 1. Particles are strongly 
forward scattering only if they are perfect, isolated, clear 
spheres [McGuire and Hapke, 1995]. Rough-surfaced, irregu- 
lar, absorbing particles have scattering patterns that are much 
less asymmetric. Hence is probably also close to 1. 

These values will be used to calculate the solid-state green- 
house effect expected for a similar material on the surface of 
an icy satellite of Jupiter. Substituting them into equation (3) 
and taking p 0 = 1 to maximize the greenhouse effect, the 
coefficient of JJ&E r k is found to have the value 3.19. Fol- 
lowing Brown and Matson [1987], it will be assumed that k — 
0.001 W nT 1 K 1 . At D = 5 AU, J v = 1360/25 - 54.4 W m -2 , 
so that T r = (irJJauY' 4 = 234 K. From equation (65) of 
paper 1, (p* Ts = p 0 (l — A h )/2tt = 0.0271, and T s = 
T r [( 1 + y* T )<p* T Jy* T Y /4 = H3 K. Also,£ r = E v - 3<J>/2</; 
it will be assumed that <f> = 0.30, but d will be allowed to vary, 
so that E v = 0.450/d m” 1 . 

When these values are inserted into (3), the temperature 
deep within the snow pack is found to be given by 


where all quantities are defined in paper 1. 

For realistic optical properties of a strongly forward scatter- 
ing, high-albedo material in the visible, the data for dry Ant- 
arctic snow measured by Schwerdtfeger and Weller [1977] will be 
used. They found that between the surface and a depth of 100 
cm the mean size of the snow particles increased from about 
0.2 mm at the surface to 1 mm at 50 cm. The intensity of visible 
light decreased according to the relation I — / 0 exp ( — x z )i 
where z is the depth below the surface of the snow and xd — 
0.03, where d is the grain size. The bulk density was p = 0.30. 
The visual albedo was 0.83, and the power of incident sunlight 
per unit area on the snow was 900 cal cm -2 d *. Since the 
irradiance of sunlight is 2800 cal cnY 2 d _1 , the cosine of the 


T(oo) = T s + 3.86 X 10 5 d - 1.26 X 10 “Vr 4 (°°) K, (4) 

where T s = 113 K is the surface temperature and d is the 
particle size in meters. This expression may be readily solved 
for T( oo). The result is plotted as the solid line in Figure 1, 
which shows the rise in temperature A T ~ T(°°) - T s versus 
regolith particle size for 10 pcm < < 1 cm. For d of the order 

of a few tens of micrometers, similar to lunar regolith, AT — 
10 K. For (7 — 1 mm, similar to the Antartic dry snow, AT — 
100 K. For (7 — 1 cm, such as in a wet, terrestrial snowfall, 
AT — 120 K. The temperature rise saturates at AT — 120 K 
for larger d. As d increases, the temperature increases, but the 
radiative conductivity also increases and limits the rise. 

These are the greenhouse temperature rises on a nonrotat- 



HAPKE: ENERGY TRANSFER MODEL APPLICATIONS 


16,835 


ing body. However, real bodies rotate with respect to the Sun, 
so that the greenhouse rise in the deep interior is proportional 
to the value of the time-dependent visible irradiance averaged 
over time. To a first approximation, this radiance is propor- 
tional to cos i averaged over the period of rotation, (cos i) ~ 
1/t t. Hence the temperature rise on a rotating body will be 
reduced by roughly a factor of 1 hr, which is indicated by the 
dashed line in Figure 1. This temperature increase does not 
exceed 40 K. 

Moreover, the low thermal inertias and inferred low thermal 
conductivities [j Brown and Matson, 1987] imply that the regolith 
particle sizes on the icy satellites are probably small, similar to 
lunar soil, so that the solid-state greenhouse effects that might 
be expected there are relatively modest, a few tens of degrees. 

Temperature rises of the order of 40 K are certainly not 
negligible and, as Matson and Brown [1989] emphasize, could 
have important implications for the interior structures of the 
bodies. However, this analysis does not confirm the potential 
increases of 400 K and more on icy Jovian satellites estimated 
by their models. Matson and Brown, of course, recognized that 
such high temperatures were physically impossible in icy rego- 
liths but argued that sublimation, melting, and the resulting 
mass transport were the processes that kept the temperatures 
low. Fanale et al. [1990, p. 203] asserted that “If (and only if) 
the effects of latent heat, mass transport and densification are 
taken into account, we can reconcile the magnitude of the 
theoretically predicted greenhouse effect with the observed 
morphology of Europa 

However, a solid-state greenhouse effect is intrinsically dif- 
ferent from an atmospheric greenhouse effect, because in the 
latter the optical depth at which visible energy is deposited is 
essentially independent of the bolometric albedo, whereas in 
the former, the two are tightly coupled. The Matson and Brown 
[1989] model did not adequately take this coupling into ac- 
count, nor did it appreciate that radiation can contribute to the 
thermal conductivity even at temperatures below 200 K. Melt- 
ing and sublimation may very well be occurring in the interior 
of Europa, as they hypothesize, but these processes are not 
necessary to limit the greenhouse temperature rise. 

3. Thermal Beaming 

The radiometric technique [Morrison and Lebofsky, 1979] is 
a major method for determining the diameters and physical 
albedos of objects that are too small to be resolved. It is 
especially important in the study of asteroids. This method 
involves comparing measured radiances in the visible and ther- 
mal IR received from the object at small phase angles. Early 
models, which assumed that surfaces radiated in the IR in 
accordance with Lambert’s law, yielded diameters that were 
systematically too large and albedos that were too small. It is 
known that many surfaces, notably the Moon’s, are not Lam- 
bertian emitters; hence an empirical parameter to correct for 
this, called the thermal beaming factor, was introduced. Most 
explanations of the nature of the beaming factor, particularly 
that of Spencer [1990], have concentrated on the effects of 
unresolved surface roughness. However, neither the visible 
reflectance nor the thermal emission from the surface of a 
particulate medium are Lambertian, so that it is of interest to 
enquire as to the importance of their contributions to thermal 
beaming. 

To estimate the magnitude of these effects it will be sufficient 
to consider only the so-called standard model, in which it is 



Figure 1. Increase in temperature due to the solid-state 
greenhouse effect, calculated from equation (4), plotted 
against regolith particle size. The solid line shows AT on a 
nonrotating body, and the dashed line shows AT on a rotating 
body. 

assumed that the body does not rotate and the regolith is in 
radiative equilibrium. The visible power scattered by a spher- 
ical body into zero phase angle is, by definition, 

P v = J v tt 31 2 ~ A pv = J t 9k 2 A, )v = J v m. 2 AJZ„ (5) 

where 91 is the radius of the body, A pv is the visual physical 
albedo, A sv is the visual spherical (Bond) albedo, and £„ = 
A p JA sv is the visual phase integral. 

The thermal power scattered by the body into zero phase 
angle is 

Pt = -M* ~ e )p ddl 

J d. 

C Till 0 - 

= s d (p) ~ T e (p)p2'n , 3i 2 sin e de, 

J o 

where si is the illuminated area on the surface of the body and 
T e is the effective temperature at which the surface radiates to 
space. On any illuminated unit area of the surface in equilib- 
rium, 

e„ — Tt( p 0 ) = hP o[ 1 - A hv ( p 0 ) 3 , (6) 

TT 

so that at zero phase where p 0 = p, 

P T - 2'n < 3l z J l , f ^^-[1 - A hv (p)]p 2 dp. (7) 

£/i 

J o 

Combining (5) and (7) gives 

^ = 2ir L ( ' — [1 - dtL. (8) 

P v A sv e H 
J o 

If the surface were perfectly Lambertian, the following 
would be true: (1) s d (p) is independent of e and equal to s h ; 
(2) A hv (p 0 ) is independent of i and equal to A xv ; (3) the 
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e (degrees) 

Figure 2. Measured effective power per unit area emitted by 
the lunar surface as a function of the angle of emergence e 
compared with theoretical models of a regolith with a macro- 
scopically smooth surface. The top points and curve corre- 
spond to emission from the subsolar point at different e. The 
bottom two curves correspond to a scan across the full Moon. 

surface radiates according to Lambert’s law, so that T, 4 (/x 0 ) == 
T*(i = 0 )/x 0 . Then the integral in (7) and ( 8 ) can be evaluated 
to give 

P t =2tt^ 2 J v (1-AJ/3, 

and 

P T /P V = (2ir/3)£„(l - AJ/A SV . 

To correct for non-Lambertian effects, equation ( 6 ) is mod- 
ified by multiplying the left-hand side by an empirical thermal 
beaming parameter tj, while keeping = s h andA hv (p 0 ) 

= A sv ; thus 

Ve»^7l(n„)=Jji„(l-AJ. (9) 

Then (7) and ( 8 ) become 

P r =2ir9l 2 J v (l ~AJ/3 V , (10) 

and 

P T /P V = (27r/3)^(l - A sv )/rp4 sv . (11) 

A value is assumed for £ v (often the lunar value), and, assum- 
ing rj is known, equation (11) can be solved for A sv and A pv — 
ZvAsv, then Sft can be found from (5). 

Requiring that (7) and (10) for P T be equal gives the fol- 
lowing expression for rj: 

1 1 ~A SV 

v = -~ . ( 12 ) 

3 f ' n a ( x\ * a 
—-—[ l-A h M\ix 2 dix 

J 0 

The value of tj based on observations of the Moon and the 
Galilean satellites is tj =* 0.75 [Lebofsky et al., 1986]. 

In order to evaluate the contribution of the intrinsic non- 
Lambertian behavior of e rf (/x) and A hv (fi 0 ) of a particulate 
medium, the Moon will be taken to be a representative body. 
Equation (69) of paper 1 was used to numerically calculate the 
radiance I T (i = 0 , e ) emitted by the subsolar point as it moves 


across the Moon’s surface during a lunation. The same equa- 
tion was used to calculate the distribution of radiance I(i = e, 
e ) across the full Moon. The equations derived in section 4 of 
paper 1 were used to calculate the various quantities in ( 12 ) 
and 17 was evaluated. The bolometric albedo of the Moon was 
taken to be A sv — 0.12 [Lane and Irvine, 1973], and the 
hemispherical emissivity e h = 0.89 [Saarietal., 1972]. A value 
was assumed for j3„, and w v was found from A sv , following 
which A h (iL 0 ) was calculated. Similarly, a value was assumed 
for /3 r , w T was found from s h , and then s d (p) was calculated. 
I T (i, e^/was then found from^,,(/x 0 ) and e d (jx). 

The calculated predictions were compared with the lunar 
thermal data reported by Sinton [1961] and Saari and Shorthill 
[1972]. The results are shown in Figure 2. The triangles are 
Sinton’s measurements of “planetary heat,” which is the power 
per unit area emitted by the subsolar point as it moves across 
the lunar surface during a lunation inferred from thermal IR 
measurements, assuming that the Moon radiates in accordance 
with Lambert’s law. However, if the Moon radiated like a 
perfect Lambert surface, the emitted power per unit area 
would be independent of e. The observed limb darkening is an 
indication of departures from Lambertian behavior. The 
model curve exhibits some limb darkening, but not nearly as 
much as the Moon. 

The observed [Saari and Shorthill, 1972] power distribution 
across the full moon is approximately proportional to ju 2/3 , 
whereas for a Lambert surface the power would be propor- 
tional to /a 1 . Figure 2 shows the power per unit area calculated 
from the model and according to /x 2/3 . The amount of limb 
darkening predicted by the model is too large. 

The thermal beaming factor predicted by the model was 
calculated by numerical integration of ( 12 ) and found to be 
tj = 0.95. 

The values of the asymmetry parameters used in calculating 
the curves of Figure 2 are /3„ = j3 r = 0. Other values were 
tried, but the resulting curves were found to be quite insensi- 
tive to this parameter. The value of A sv was varied between 
0.09 and 0.15 and of e h from 0.85 to 0.95, again with only minor 
effects. 

Hence although the angular emitting and scattering proper- 
ties of particulate media are not Lambertian, the departures 
are insufficient to explain the observations. In particular, they 
account for only 20% of the departure of rj from unity. The 
Moon exhibits a surge in visible brightness at full Moon, and it 
may be speculated that a similar effect might occur in the IR. 
However, the visible opposition effects from planetary rego- 
liths result from either or both coherent (coherent backscatter) 
or incoherent (shadow hiding) phenomenon in which the di- 
rection of the incident radiation plays a special role. Since 
there is no external source of thermal radiation, analogous 
phenomena cannot occur in the IR. Thus the other 80% of tj 
is most likely due to macroscopic roughness [Spencer, 1990]. 

4. Spectral Emissivity 

The location and contrast in features of emissivity spectra 
are important for infrared remote sensing and have been the 
subject of considerable recent attention. Conel [1969] and Lo- 
gan and Hunt [1970] identified maxima in the emissivity rest- 
strahlen spectra of particulate materials with the Christiansen 
wavelength and pointed out that this wavelength was diagnos- 
tic of composition. Since then these maxima have been inves- 
tigated in detail by several others, including Salisbury [1993], 
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Salisbury and Walter [1989], B. Henderson and B. Jakosky 
Near-surface thermal gradients and mid-IR spectra: A new 
model including scattering and application to real data, (sub- 
mitted to Journal of Geophysical Research , 1995), and Wald and 
Salisbury [1995]. A particular concern has been the effect of 
subsurface temperature gradients on the wavelength of the 
maxima. Henderson and Jakosky [1994, also submitted manu- 
script, 1995] concluded that such temperature gradients had 
only minor effects. 

The formalism developed in previous sections of this paper 
will be used to investigate this question further. Other authors 
used the spectra of real minerals. However, such spectra usu- 
ally consist of multiple overlapping bands, which tends to ob- 
scure important properties of the spectra. Hence the Lorentz 
model of absorption bands will be used to generate theoretical 
emissivity spectra of media in radiative equilibrium with sun- 
light in order to clarify and elucidate the nature of the features 
of a strong reststrahlen band seen in emissivity. Lorentz theory 
[e.g., Wooten, 1972] is widely accepted as providing a reliable 
description of such bands [e.g., Spitzer and Kleinman, 1961]. 

According to the Lorentz model, the complex dielectric con- 
stant (also called the dielectric function) in the vicinity of a 
strong absorption band is given by 

K e = K r +iK„ (13a) 


the following expression is a good approximation for w x when 
ray theory is applicable, 

= ( 15 ) 

where 

a — 4 tt nJX (16) 

is the absorption coefficient and S e is the coefficient for reflec- 
tion of radiation externally incident on the surface from all 
directions, 


S e = [/?,(#) + R 2 (d)] cos d sin ■& d-d, (17) 


where /?,(#) and R 2 ($) are the Fresnel reflection coefficients 
[Born and Wolf, 1980], 


[cos d — G]] 2 + G 2 
R 'W = [cos d + G,] 2 + GV 


(18a) 


[(nj — nj) cos d — Gj] 2 + [2 n,n, cos d G 2 ]“ 
Rii^) [(nj — nj) cos d + G,] 2 + [2 ngt, cos d + G 2 ] 2 

(18b) 


where 


„ , vlivl-v*) 

Kr ~ K+ (vl- v 2 ) 2 + H V ’ 


K,= 


(vl-v 2 ) 2 + SV’ 


(13b) 

(13c) 


G,(d) = {\[((nj - nj - sin 2 d) 2 + 4 njnj)' 12 
+ (nj - nj - sin 2 d)]} 1/2 , 

G 2 (d) - {\[((nj - nj - sin 2 d) 2 + 4 njnj) 112 
- (nj - nj - sin 2 d)]} 1/2 . 


(18c) 


(18d) 


i = V— T, n c is the continuum real part of the refractive index 
outside of the band and v is the frequency; the band is char- 
acterized by three parameters: the frequency of the band cen- 
ter v Q , the plasma frequency v p , and the collision frequency a. 

By definition, the dielectric constant is the square of the 
complex refractive index, 

m—n r +in h (14a) 

so that 

K c = K, + iKj = nj. - nj + i2n r n„ 
which may be solved for n r and n i to give 


n,= [\m + Kr)' n + K ,]}'« 

(14b) 

«, = !jiW + rf)' n - K,]) vl . 

(14c) 


In order to calculate the emissivity, a model for the wave- 
length-dependent single scattering albedo w A and angular scat- 
tering function p x (g) of the medium is needed. It will be 
assumed that the particles are all equant and of a single size d 
large compared with the wavelength A. In that case, Q EX ~ 1 
and the single scattering albedo w x is equal to the particle 
scattering efficiency Q sX . Mie theory will give incorrect results 
if it is used to calculate w A because it applies only to an isolated 
sphere and includes diffraction, which is inappropriate for par- 
ticles in a close packed medium [Hapke, 1993]. Since the par- 
ticles are large compared to the wavelength, geometrical op- 
tics, or ray theory, is appropriate. Hapke [1993] showed that 


In (15) the contributions of different phenomena to particle 
scattering are explicitly separated. The first term of (15) de- 
scribes radiation that has been reflected by the surface of the 
particle, while the second term describes light that has been 
refracted into the particle and partially absorbed there or scat- 
tered back out from within its volume. All orders of internal 
scattering are included. 

Although transparent spherical particles are strongly for- 
ward scattering, highly absorbing ones are not. Furthermore, 
McGuire and Hapke [1995] have shown that any departures 
from a perfect sphere, such as irregular shape or surface 
roughness causes the particle phase function to be more iso- 
tropic. Hence, for simplicity, it will be assumed that the parti- 
cles scatter isotropically in both the visible and IR, so that the 
particle phase function p x (g) = L 

Since the particles are large compared to the wavelength, the 
extinction coefficient E x will not vary appreciably in the ther- 
mal IR so that t a =* t t . Then for isotropic particles the 
thermal radiance emerging from the surface of the medium at 
wavelength A and angle e is [Hapke, 1993] 




w x <p x (r T ) + — U X (T) 


e-^dTj/p, (19) 


where <p A (r r ) is the average isotropic spectral radiance per 
unit wavelength interval, 

7a = - w A > (20) 

U X (T) - (2Trh l)C{) /\ 5 )(e ,nMT - D' 1 ( 21 ) 



16,838 


HAPKE: ENERGY TRANSFER MODEL APPLICATIONS 




Figure 3. Lorentz model spectrum of a reststrahlen band 
with parameters v p /v 0 = 1.00, u/v Q = 0.01, n c = 1.50. The 
band center is at A 0 : (a) dielectric constant and (b) refractive 
index. 


A 0 = cjv 0 = 10 (xm, d! A 0 = 10, w„ = 0.36, j3„ = p T = 
0. The spectra of the resulting dielectric constant and refrac- 
tive index are shown in Figures 3a and 3b, respectively, and the 
single scattering albedo is shown in Figure 4. 

Four cases of interest were studied: q = 0.26, /x 0 = 1, 
T r — 524 K, w T - 0.52, corresponding to the subsolar point 
on the Moon viewed normally; q - 0.26, jx 0 = cos 60°, T r = 
524 K, w T = 0.52, corresponding to a point on the lunar 
surface 60° away from the subsolar point and viewed normally; 
q = 1.067, p o = 1, T r - 839 K, w T = 0.52, corresponding 
to the subsolar point on Mercury viewed normally; and q = 
0.023, p, 0 = 1, T r = 234 K, w T - 0.49, corresponding to 
the subsolar point on a Jovian satellite viewed normally. These 
parameters were chosen in order to vary the subsurface tem- 
perature gradient over a large range, as shown in Figures 3-7 
of paper 1. This gradient is largest for the subsolar point on 
Mercury, small for a satellite of Jupiter and negligible for i = 
60°. The values of w T were found by weighting w A by a Planck 
function with T = T r and integrating over wavelength, as 
discussed in paper 1. 

According to equation (71) of paper 1, the medium may be 
treated approximately as emitting like a surface at the temper- 
ature T s . Hence the emitted spectral radiance / A (/x) was con- 
verted to an apparent spectral directional emissivity s dK ([x) by 
dividing by U K (T S ) hr. The resulting apparent spectral emis- 
sivities are shown in Figure 5. Also indicated in Figures 3-5 is 
the location of the Christiansen wavelength. Figures 3-5 illus- 
trate several interesting conclusions. 

First, the emissivity minimum is not at the true center of the 
reststrahlen band, A = A 0 , but is displaced well toward shorter 
wavelengths. The reason can be seen by inspection of Figures 
3-5. The emissivity minimum corresponds to the maximum in 
w A , which occurs where the real part of the refractive index n r 
« 1, but the imaginary part n, is still relatively large, so that 
the volume scattered term of w k is negligible. To a rough 
approximation, 

(n r ~ l) 2 + nj 

c ( n r + l) 2 + nf ’ 


is the Planck function, h 0 is Planck’s constant, c 0 is the speed 
of light in vacuum, and k 0 is Boltzmann’s constant. 

The wavelength-dependent boundary layer approximation 
for the radiance is given by equation (75) of paper 1, which for 
isotropic scatterers (£ A — 1) is 


<p A (r T ) - -j- U k {T) 


9a£ 

7a 


~2yjJT 


where 


<Pas “ 


7a 


1 


U k (T s ), 


1 + 7a 7T 

T(r r ) - T r T*(r* T ), 


( 22 ) 

(23) 

(24) 


T*(t* t ) is given by equation (60) of paper 1, which describes 
the temperature regime in a medium in radiative equilibrium 
with sunlight, and 


T,T\ 


(25) 


where T* is given by equations (61) and (56) of paper 1. 

This system of equations was converted to wavelength as the 
independent variable and evaluated numerically with the fol- 
lowing parameters: n c = 1.50, v Jv 0 = 1.00, H/n 0 = 0.01, 



A/A 


Figure 4. Single scattering albedo spectrum of a spherical 
particle with d — 10A o and the refractive index shown in 
Figure 3b. Solid line is single scattering albedo; dashed line is 
contribution of the surface-reflected component, S e . 
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hence, if n r « 1, S e is close to 1, w A has a maximum, and 
e dx (p) has a minimum. 

Second, the short-wavelength maximum in apparent emis- 
sivity is close to, but not exactly at, the Christiansen wavelength 
A c , which, by definition, occurs where n r - 1.00. This is 
clearly seen in Figure 6, in which the continuum refractive 
index n c is varied in order to change A c . When n c is large, the 
emissivity maximum is very close to A c . However, as n c de- 
creases, A c is increasingly displaced toward shorter wave- 
lengths away from the emissivity maximum. 

The traditional argument [Conel, 1969; Salisbury, 1993] is 
that the medium is transparent there because particle surface 
scattering is negligible; hence radiation from deeper, hotter 
layers of the medium reaches the surface, increasing the ap- 
parent emissivity. However, Figures 4 and 5 show that the 
maximum in the emissivity actually occurs at the place that 
Hapke [1993] has called the first transition minimum in w A and 
in reflectance. As A increases near this minimum, n i also in- 
creases rapidly, until the particle is nearly opaque and the 
single scattering albedo changes from being dominated by vol- 
ume-scattered light to being dominated by strong surface- 
scattered light. Even though the surface scattering is small, the 
medium is not transparent because the absorption coefficient is 
not negligible there: at A = A c , ad = 1.26. 

Fortuitously, this first transition minimum always occurs 
close to A c . Hence it is appropriate to continue to refer to the 
associated emissivity maximum as the “Christiansen maxi- 
mum,” especially since this terminology is now ingrained in the 
literature. However, this maximum would occur even if n r did 
not go through 1.00 anywhere and even in the absence of any 
subsurface temperature gradients. 

Third, there is another maximum on the long-wavelength 
side of the emissivity minimum. Salisbury and Walter [1989] 
have called this the “transparency feature” and attributed it to 
a second region of relative transparency of the medium be- 
tween two reststrahlen bands. However, Figures 4 and 5 show 
that this maximum is actually caused by a second transition 
minimum in vt> A , where the particle albedo changes from being 
dominated by strong surface scattering to being dominated by 
volume scattering as n t decreases with increasing wavelength 



Figure 5. Apparent emissivity spectra of a regolith com- 
posed of isotrotropically scattering particles with the spectral 
single scattering albedo shown in Figure 4, in equilibrium with 
sunlight, and illuminated and observed under different condi- 
tions. 


n =1.3 C 1.5 L7 



Figure 6. Emissivity spectra, showing the effect of varying 
the continuum real refractive index n c . The arrows indicate the 
wavelengths of the corresponding Christiansen wavelengths 
A c . Note the displacement of A c from the emissivity maxi- 
mum. 


[Salisbury, 1993; Hapke, 1993]. It occurs even in the absence of 
a second absorption band. The wavelength at which the trans- 
parency maximum occurs is much less sensitive to n c than the 
Christiansen maximum. 

Salisbury and Walter [1989] have noted a strong correlation 
between the wavelengths of the emissivity minimum and the 
two maxima. Such a correlation is to be expected because all 
three are associated with the reststrahlen band. Thus all should 
be reliable indicators of composition, especially if they can be 
resolved in the same spectrum. The widths of the two maxima 
should depend on particle size: as the particle size increases, 
both the Christiansen maximum and the transparency feature 
should increase in width. 

Fourth, a large subsurface temperature gradient causes the 
Christiansen maximum to shift toward shorter wavelengths 
because of the increasing radiance at short wavelengths from 
the Planck function. However, the amount of the shift is small. 
This confirms the previous conclusion of B. Henderson and B. 
Jakosky (submitted manuscript, 1995), and is also consistent 
with laboratory measurements by Salisbury et al. [1994b], The 
reason for the lack of sensitivity is that the particles are not 
particularly transparent there, as previously thought. 

Fifth, the main effect of the subsurface gradient is to raise or 
lower the emissivity nearly uniformly over the entire wave- 
length region. This is why the apparent emissivity exceeds unity 
for some of the curves. However, in practice, this should have 
little effect. Because of calibration difficulties, absolute emis- 
sivity is seldom measured. Usually, emissivity is arbitrarily nor- 
malized to unity at the Christiansen maximum. 

Also shown in Figure 5 is the effective spectral directional 
emissivity, 

^ £a7a [1 + 2£aF-] + /3a[1 ~ 2t a /£a ~ 

“ 1 + • 1 - J3a 

(26) 

which can readily be derived by obvious wavelength-dependent 
extension of equation (71b) of paper 1. Equation (71b) of 
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paper 1 was derived using the boundary layer approximation, 
and also assuming a constant subsurface temperature. In Fig- 
ure 5 it was assumed that /3 X = 0, in which case, equation (26) 
is 

1 + 2/x 

• ( 27 ) 

Note that its shape is virtually identical to those of the other 
curves in Figure 5. Hence, except in the unusual circumstance 
that absolute emissivity is required, the approximate expres- 
sions (26) or (27) may be reliably used to infer or calculate the 
IR parameters of the medium. 

5. Conclusions 

Using a new model [ Hapke , this issue] that simultaneously 
solves the equation of radiative transfer for visible and thermal 
radiation and the heat equation, three problems of interest in 
the thermal properties of planetary regoliths were studied. 

The solid-state greenhouse effect is shown to be self-limiting 
because of the strong correlation between optical absorption 
and albedo and because radiative thermal conductivity is im- 
portant in a porous regolith, even at low temperatures. If the 
optical properties of terrestrial snow are representative of 
those of the regoliths of satellites in the outer solar system, the 
maximum rise in temperature that may be expected due to a 
greenhouse effect is about 40 K. 

The contribution to thermal beaming caused by non- 
Lambertian albedos and emissivities of a macroscopically 
smooth regolith was estimated to be able to account for a 
thermal beaming parameter of about 0.95, whereas the ob- 
served value is about 0.75. Since coherent backscatter or 
shadow hiding opposition effects cannot occur in thermal ra- 
diation from smooth regoliths of small particles, it is likely that 
the remainder of the beaming factor is caused by macroscopic 
roughness. 

The nature of features in the emissivity spectrum of a rest- 
strahlen band was investigated using an artificial absorption 
band generated by the Lorentz equation. It was found that the 
peak that occurs on the short wavelength side of a spectral 
emissivity band has been incorrectly identified as denoting the 
Christiansen wavelength. This was a natural error, because the 
wavelength of the maximum is close to, but not exactly at, the 
Christiansen wavelength. The emissivity peak corresponds to 
one of two transition minima in reflectance, where the radia- 
tion scattered by a soil particle changes from volume scattering 
to surface scattering. The other transition minimum occurs on 
the long wavelength side of the band and corresponds to the 
transparency feature seen in emissivity. Because the particles 
of the regolith are not transparent at either peak, subsurface 
temperature gradients have only a minor effect on the wave- 
lengths of their locations. Hence the wavelengths of the two 
peaks and the included minimum should be reliable diagnostic 
indicators of composition. 
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